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• Titre en ffancais: Relation de dispersion d'une couche de ferrofluide de viscosite et 
d'epaisseur quelconques sous champ magnetique normal; Regimes asymptotiques. 



• Classification Physics Abstracts: 

1. 47.20. Hydrodynamic stability and instability. 

2. 75.50.Mm. Magnetic liquids. 

3. 47.35. Waves. 

• Resume en Francais: 

Nous avons calcule la relation de dispersion des ondes de surface dans une couche de 
ferrofluide d'epaisseur et de viscosite quelconques, soumis a un champ magnetique 
normal a sa surface (instabilite de pics de Rosensweig) . Cette relation montre que 
le champ magnetique critique et le vecteur d'onde critique de 1'instabilite dependent 
de Pepaisseur de la couche de fluide. La relation de dispersion a ete simplified pour 
quatre regimes asymptotiques: couche epaisse ou mince et comportement visqueux 
ou inertiel. Nous avons calcule les valeurs critiques de 1'instabilite dans ces quatre 
cas. Nous montrons qu'un parametre typique du ferrofluide permet de savoir dans 
quel regime, visqueux ou inertiel, se situe le ferrofluide pres du seuil de 1'instabilite. 
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Abstract 

We have calculated the general dispersion relationship for surface 
waves on a ferrofluid layer of any thickness and viscosity, under the in- 
fluence of a uniform vertical magnetic field. The amplification of these 
waves can induce an instability called peaks instability (Rosensweig 
instability). The expression of the dispersion relationship requires that 
the critical magnetic field and the critical wavenumber of the insta- 
bility depend on the thickness of the ferrofluid layer. The dispersion 
relationship has been simplified into four asymptotic regimes: thick 
or thin layer and viscous or inertial behaviour. The corresponding 
critical values are presented. We show that a typical parameter of the 
ferrofluid enables one to know in which regime, viscous or inertial, the 
ferrofluid will be near the onset of instability. 



1 Introduction 

It is known since the experiment performed by Cowley and Rosensweig (ffl, 0) that a 
normal magnetic field has a destabilizing influence on a flat interface between a magneti- 
zable fluid and a non magnetic one. Above the magnetic induction threshold H cr u, the 
initially flat interface exhibits a stationary hexagonal pattern of peaks (Rosensweig insta- 
bility). The hexagonal pattern becomes a square pattern above another critical magnetic 
induction H' crit > H crit j|. 

In this paper, we want to obtain a general dispersion equation (ferrofluid layer of any 
thickness and any viscosity under the influence of a uniform vertical magnetic field) by lin- 
earizing the equations governing the problem. This allows to calculate the critical values 
of the instability: wavenumber and magnetic induction. Nevertheless, linearized equa- 
tions are inadequate for the full description of either phenomenon: the symmetry of the 
developed wave array (hexagonal or square) and the wave amplitudes may be determined 
only from nonlinear equations. It has been treated by Gailitis in 1977 ||. A generalized 
Swift-Hohenberg equation constitutes a minimal model to account for the formation of 
hexagons as well as squares. In such an envelope equation, for which the quadratic term 
is sufficiently small, hexagons are to be expected from the modeling B. Furthermore, 
a weakly nonlinear analysis involves that the dynamics of the patterns depends on their 
symmetry B. 



We shall remember previous calculations of dispersion equations in section (1.1) and sec- 
tion fll.2p will put in evidence the analogy with the electric case. Indeed, a layer of liquid 
metal under a normal electric field develops a similar instability of peaks (Taylor cones 
H). We applied to the magnetic case the previous analyse of Hynes and Limat || 
for the gravitational amplification of capillary waves (Rayleigh- Taylor instability) and of 
Neron de Surgy et al. || for the electric amplification of capillary waves (Taylor cones) 
H . The general dispersion relationship that we obtain in section (||) leads to the fact that 
the magnetic induction threshold and the critical wavenumber depend on the thickness 
of the ferrofluid layer fici|| . We derive the asymptotic behaviour in various regimes and 
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give analytical dispersion equations in the case of thick or thin, inertial or viscous layers. 
From a general approach, we then recover three previously known results (thick-inertial 
D, thick -viscous [jnj, thin-viscous obtained from various approaches. Furthermore, 
we give an explicit expression of the critical magnetic induction in the case of a thin layer 
of ferrofluid: it differs from the one of the thick case. We also demonstrate the influence of 
the ratio l v /l c , where l c is the capillary length and l v the viscous length of the ferrofluid: 
depending on the value of this ratio, the ferrofluid will have a viscous or inertial behaviour 
near the onset of the instability. 



1.1 Previous approaches 

Ferrofluids or magnetic liquids are permanent, colloidal suspensions of ferromagnetic par- 
ticles (100 A) in various carrier solvents. Fluid instabilities can arise with these liquids, 
especially surface instabilities. In the case of the peaks instability, the study of the linear 
stability of the interface, by Cowley and Rosensweig Q, between an half-infinite and in- 
viscid ferrofluid of density p and the vacuum, with normal modes of perturbation, leads 
to the following dispersion relation: 

1 + Mo/M 

where s is the growth rate of perturbation, k (the modulus of k) is the horizontal wavenum- 
ber, /i is the magnetic permeability of the ferrofluid, 7 is the interfacial tension ferrofluid- 
vacuum, Hq the modulus of the normally applied magnetic induction and g the gravita- 
tional field. 

The critical values of the instability are then: 



J Hc ™t - (/*/«>- W^7 (2) 
I k crit = y/pg/'f = k c = l/l c 

where k c is the capillary wavenumber. 

Above the threshold of linear stability, a broad band of wavenumber is theoretically found 
unstable. It is seen from other instabilities that the mechanisms of wavenumber selection 
are very complex, including the effects of sidewalls, time-dependent effects in cellular 
structures, axisymmetric stucture in the case of Eutectic solidification, the nonlinearities 



of the problem, etc.. 1 13 . A common choice has been to consider then that one wavenumber 
k m is selected. This wavenumber corresponds to the maximum growth rate s m of the small 
undulations, solution to the following equation dks\u=const = 0. This prediction has only 
be checked when Ho increases very quickly from a subcritical value to H C rit- In most 
cases, the wavenumber remains equal to k cr n above the onset, whatever the symmetry 
is or the experimental procedure used to increase the field (continuous increase jj], field 
jumps (U or alternating field (l5) at different frequencies). This encouraged D. Salin, in 
1993, to take into account the effect of the ferrofluid viscosity, adapting the Landau and 
Lifshitz approach Jig] . His result leads to a wavenumber k m equal to k cr i t , whatever the 
field is |ilf| . Simultaneously, Neron de Surgy et al. M obtain the same result for metal 
liquids. As seen above, the experimentally selected wavenumber is not obviously k m nor 
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kcrit p7| - We shall see in section (4.1.3) that the description of the unstable behaviour of 
the ferrofluid (viscous or not) near the onset can be determined by the value of the ratio 
l v /l c . When this ratio approaches zero, an inertial description is actually sufficient. 



1.2 Analogy with electro-capillary instability 

In 1993, G. Neron de Surgy et al. studied the linear growth of electro-capillary instabilities 
in the very general case where the viscosity of the fluid and its thickness are of any value 
H , following the studies of Hynes and Limat || . The study of this instability is similar 
to our study: an electric field that is applied normally to the free surface of a conducting 
fluid (mercury, for example) has a destabilizing effect on this interface. In their paper, 
they derived the asymptotic behaviour in various regimes and gave dispersion equations 
in the case of thin or thick, inviscid or viscous films: we present a similar approach for 
the magnetic case (but note that substituting E for H in the electric dispersion equations 
does not result in the correct magnetic dispersion equations). Furthermore, we note that 
the critical electric field and the critical wavenumber of the electro-capillary instability 
does not depend on the thickness: While E = in the whole layer of mercury, there is 
penetration of the magnetic field in the ferrofluid and the distribution of this field obviously 
depends on the layer thickness (boundary conditions). 

2 Characteristic scales for various asymptotic 
regimes 

We consider an incompressible and viscous film of magnetic fluid of density p, dynamic 
viscosity rj, kinematic viscosity v and magnetic permeability /i, that occupies the space 
between z = and z = —a, in the vacuum. The geometry is supposed to be infinite for 
both x and y and the medium above and below the ferrofluid is also infinite, so that we 
consider a two dimensional (x, z) problem (Figure |^). 

If we want to introduce the thickness effects, we have to consider a horizontal length scale 
lx — A/2-7T (with A = 27r/fc), and a vertical length scale l z . The value of l z will be taken as 
the lowest value between the film thickness a and X/2tt. 

For the study of the viscosity effects, we shall introduce the Reynolds number (Re = inertia 
forces/viscous forces). The Reynolds number is: 

If! 

" v[{l/l x ) 2 + (l/l z ) 2 } 

where v is the fluid velocity and A the Laplacian operator, 
-for a thick film, where the vertical scale is X/2n, Re = -Jtt. 
-for a thin film, where the vertical scale is a, Re = ^ a . 

If Re ^> f, we may neglect viscosity and call the film inertial, and if Re -C 1, we may 
neglect inertia and call the film viscous. 

This study results in four asymptotic behaviours (Table Q, 0, ||, ||). 



Re 



pd t v 



rjAv 
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3 Equations of the problem 



We shall present the equations governing the system shown in Figure (^). Linearity al- 
lows us to treat the different Fourier components separately. We shall denote £(x,z) the 
deformation of the interface, v(x, z) the velocity of the fluid, H = Ho + h the perturbed 
magnetic induction and n the unit vector normal to the interface where n ~ (d x S,, — l)u z 
at first order. We consider p(H) = fj,(H cr ii) because our analysis is valid near the onset 
(see section (||)). 

Local equations governing the motion of the magnetic fluid are: 

V • v = (continuity equation) 

p[d t v + (v • V)v] = —V p + i]Av + pg (Navier-Stokes equation) 

Other governing relationships are the Maxwell equations: 

V x H = (no charges, no currents) 

V • H = 

The boundary conditions (where [[X]] = value of X above the interface — value of X 
under the interface) are given by: 

dt£, = v z — v x d x £, — Vydyt; at z = £, (free surface condition) 

— [[p]]ni + [[To; + a' ik ]]nk — {j/R)rii = at z = £, (sress balance at the interface) 

[[n • /xH]] =0 at z = £ and z = —a 

[[n x H]] =0 at z = £ and z = —a 

v z — at z = —a 

v x = at z = —a 

where is the stress tensor given by (0], [§]): 

Tik = pHiHk — ^H 2 Sik 
and a' ik is the viscous rate-of-strain tensor given by: 

v'ik = v(d Xk Vi +d Xi v k ) 
and is the curvature of the interface (positive if directed towards the fluid): 

d 2 d 2 
^ dx 2 dy 2 ^ 

We seek solutions with the following form: 

A = sRe[A(z)exp(st - ikx)] 

for £, h Xl h z , v x , v z . At first order, we get a system of nine linear equations in nine 
unknown. In order to find a non trivial solution, the following determinant has to equal 
zero. Similar calculations can be read in the electro-capillary case H. 
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with q 2 = k 2 + sjv. 

We shall now use dimensionless values with capillary scale as reference as seen in Table 
||. We introduce the parameter / = (l v /l c ) 3 ^ 2 (denoted d in g) where the viscous length 
U = u^g- 1 / 3 (§,§). We obtain that q 2 = k 2 + f 

We also denote $ = HQ/(H^ k ) 2 where = H crlt (section (|j~Tl)). The dispersion 



relation s = s(Hq, k) is an implicit equation and its dimensionless form can be written: 

Aak 2 (k 2 -\~ Q 2 } 

4qk 3 (q- fc coth(ga) coth(fca)) - (fc 2 + <? 2 ) 2 (gcoth(fca) coth(qa) - k) 1 ' '" 



sh(qa)sh(fca) 



= _L(fc + fc 3 - 2$fc 2 1 + ^J^ )(gcoth(go) - fccath(fca)) (3) 

■' 1 + 1 + F(ka) 



where 

F(ka) ^ e 



(1 + aV^o coth(fca))sh(fea) 
We can develop equation (||) when s tends to (s = is a root of equation (||)). We 
consider then that 6k — q — k ~ tt^-t and we obtain that fc + fc 3 — 2<&fc 2 ~ dfc. 

-L+1 + F (fca) 

<5fc = (i.e. s = 0) leads to the curve of marginal stability. We obtain that k cr n and 
§ C rit depend on a (M, |Io| ). We shall now study the dispersion relationship in asymptotic 
cases. 



4 Asymptotic behaviour of the dispersion re- 
lation 

4.1 Thick film 

The equation (|3|) can be simplified in the regime a> A (fca 3> 1). 
4.1.1 Thick-inertial film 

This corresponds to the case where Re ^> 1 i.e. j^p" ^ 1 (dimensionless Reynolds). 
The equation (0) becomes the well-known Cowley and Rosensweig's result El: 
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.s 2 = -fc 3 + 2$fc 2 - fc 



(4) 



where $ cr ;t = 1 and k cr n = 1. As was shown previously, the wavenumber corresponding 
to the maximum growth of perturbations, is found to be field-dependent: 



Let us consider the parameter e so that $ = 1 + e and the parameter 6k, where k = 1 + 6k 
(6k depending on which wavenumber(s) is (are) actually selected). Near the onset, we 
develop the growth rate s(<!>, fc) into a series in e and 6k. We found at lowest order in e 
and 6k (e and 6k independent): 



As seen above, a band of wavenumbers of width e 1 / 2 near threshold is unstable |DJ. This 
means that the order of magnitude 6k 2 is at most e so we can take s 2 to be of order e. 
If fc m is actually selected, we check that 6k 2 is of order of magnitude e 2 (which is less 
than e). If the wavenumber selected is fc = 1 + \[2s (we remain on the curve of marginal 
stability), we have checked that a similar analysis is adequate, by making calculations at 
the following order. The validity conditions, at the lowest order, for this regime can be 
summed up as follows: 




s 2 ($ = 1 + e, fc = k crit + 6k) = s 2 (e, 6k) =2e- 6k 2 




a > 1 (fca » 1, thick film) 

f 2 < ^ £ (2m ^ 1 and 2ft of order of £l/2 //' inertial film) 



4.1.2 Thick-viscous film 



relation (§, @): 



In this case, Re <C 1 i.e. 




<C 1. The equation 



leads to the viscous-dominated 




(5) 



where 4> C rit = 1, fccrit = 1 and k m = k cr u = 1, whatever $ is. 

Near the onset of instability, the growth rate becomes after development: 



s(e,6k) = (2e~6k 2 



)/(2/) 



We can then take s to be of order e/f. The validity conditions become: 



{ 



a > 1 



P » e 



The condition of viscous regime are always true at the onset. 
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4.1.3 Range of the asymptotic regimes 



As seen above, the cross-over between the inertial and the viscous regime appears when 
e = f 2 where e <C 1 (Figure ||). This relation yields the following consequences: 

• Strictly at the onset of instability (e = 0), the condition of inertial regime is never 
reached: the ferrofluid has a viscous behaviour (/ is finite). 

• By increasing the magnetic induction, we reach the inertial regime above £ — f 2 (while 
keeping e <C 1). 

For example, with the ferrofluid APG 512 A (/ ~ 0.27), the inertial regime takes place 
when H ~ fljSJtt*(l + §) = H%% k {l + £) i.e. H ~ 1.04 H^ff. 

In Table pi we present different values of the parameter /. With the ferrofluid EMG 507, 
the inertial regime is reached when H ~ 1.00002 H^f : we can then consider that the 
description thick-inertial is adequate, whatever H is. 

It appears then, that the regime viscous or inertial, near the onset, depend on the type 
of ferrofluid. One also could use a ferrofluid, for example APG 067 or APG 314, that en- 
ables to remain in the viscous regime far enough from the onset. Physical data of various 
ferrofluids are presented in Table H. 



4.2 Thin film 

It seems obvious that the thick case is easier to experiment than the thin case. In order to 
stick to an opinion, the typical ferrofluid APG 512 A has a capillary length l c ~ 1.7 mm 
so that for the thin regime, the condition a <C l c is quite difficult to experiment. The thin 
regime could be easier to attain if we could increase l c and this may be achieved under 
microgravity conditions. Typical conditions of microgravity experiments lead to a value of 
10~ 2 9 (1 9 — 9.81 m/s 2 ) in parabolic flys pl| (Z c ~ 17 mm) and 10 -6 g in orbital stations 
[ fl9| (l c ~ 1.7 m). In order to make an experiment, one can choose a value of g that enables 
to have a value of l c large enough but not too much (the distance between the peaks is 
also proportional to l c ). Brand and Pettit made an experiment with a 8 mm deep layer 
on a parabolic fly 0] (As the capillary length of their ferrofluid is l c = 15.5 mm in the 
low gravity phase of the parabolic fly, they were not actually in thin nor thick regime) . 

4.2.1 Thin-inertial film 

If the Reynolds number Re ^> 1 i.e. ^j— ^> 1, we obtain the following equation)]: 

s 2 = a(-k 4 + 2$fc 3 (±±l^j _ fc 2) ( 6 ) 

Let us consider the function c(/i) = l^MlE^ where c(/i) < 1. At the onset of the instability, 
$ = 1+ 2 (> 1) and k crit = 1. When $ > 1+ 2 the wavenumber of maximum growth 
rate is given by: 

k m = ^(3$c + v^^c 2 - 8) 
1 see [^l| for the electric analogy 
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The asymptotic critical magnetic induction in the thin case is different from the one of the 
thick case and always larger. In order to quantify this difference, we have to solve implicit 
equations in section (||). 

We denote e as e — $c — 1. Near the onset of instability, the growth rate can be developed, 
at lowest order in e and 8k, and leads to: 

s 2 ($c = l + e,k = k crit + 5k) = s 2 (e, 8k) = a(2e - 8k 2 ) 

We obtain then that s 2 has the order of magnitude ae. 
The validity conditions are: 

a < 1 (thin film) 

f 2 /a 5 <C £ (inertial film) 

4.2.2 Thin-viscous film 

When the Reynolds number JJe« 1 i.e. ^j— <C 1, the equation (||) becomes [ p"2| : 

S = ^(-k i + 2^( 1 -±^)-k 2 ) (7) 

At the onset of the instability, $ = 1+ 2 ^ and k cr u = 1. When $ > 1+ 2 y^ , the 
wavenumber of maximum growth rate is given by: 

k m = i(3$c + \/9$ 2 c 2 - 8) 

Near the onset, we develop the growth rate in e and <5fc: 

s($c = 1 + £, A = k crit + 8k) = (a 3 /3/)(2e - Jfc 2 ) 

or also s of order of magnitude (a 3 //)e. 
The validity conditions are: 

a < 1 (thin film) 

/ 2 /a 5 S> £ (viscous film) 

The dispersion relations of the different cases are presented on Figure ||. 

4.2.3 Near the onset of the instability 

In the thin case, we can notice that: 

• At the onset of the instability, the thin film of ferrofluid has a viscous behaviour. 

• The cross-over to the inertial regime is reached when e — f 2 /a 5 (Figure ||). 

In the case of a thin layer, tends to 0: it results in the fact that a thin film has a viscous 
behaviour, at the onset and far above. With the ferrofluid EMG 507 (/ ~ 0.0068), we 
reach the inertial regime when H ~ + = #*£| t n (l + / 2 /(2a 5 )) i.e. H ~ 3 iJ*^ n 

(a = 0.1). 
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5 Implicit equations for the asymptotic crit- 
ical fields 

As the magnetic permeability depends on the magnetic field, each asymptotic value of the 
critical magnetic induction (thick or thin) satisfies an implicit equation. 

• In the thick case, the onset is given by equations (0): it is then necessary to know 
Hr{H = H^l^ k ) to calculate the critical values. In calculations of section (||), the function 
Hr(H) is a constant equal to fj, r (H*^ k ) (since our linear approach is valid near the onset). 
This leads to the implicit equation: = f(ii r {H t £% k )),J?jj^{pg~() 1 l A where f(fi r ) = 

^/ 1+1 J^ r m xhc function /(/i r ) i.e. H^it k varies rapidly as a function of fi r . 

• In the thin case, the implicit equation can be written: 



Tjthtn / Tjthick / q\ 

H crit - \ I -\ | 1 ^crit l°j 

V 

where \i T is the relative permeability of the ferrofluid. H^fr k is determined by the equa- 
tion (^) but note that in equation (||), H l c ^ k is a constant. 

The question arises as to when we can reach the asymptotic value y/2 of the ratio i?*^ ? ™/7J'^ t cfc 
in order to put easily in evidence experimentally the difference between the asymptotic val- 
ues of the critical magnetic inductions. The Langevin's classical theory has been adapted 
to yield the superparamagnetic magnetization relationship between the applied field Hq 
and the resultant magnetization M of the particle collection. For a colloidal ferrofluid 
composed of particles of one size, we have: 

M = M sat fcoth(a) - ^ = L(a) (9) 

whith a = ^P-, where m is the magnitude of the magnetic moment of a particule, L is 
the Langevin function, M sat the saturation moment of the ferrofluid and k the Boltzmann 
constant. When the initial permeability is appreciable, it is no longer permissible to 
neglect the interaction between the magnetic moments of the particles (J2^], |Q, |f24| ) 
and equation (^) has to be modified. We shall consider equation (^|) as a good first 
approximation (while \i is n °t much larger than unity). This leads to the expression of 
the relative magnetic permeability: 

This function of H decreases from /i r (0) = Xi + 1 to 1. With the ferrofluid APG 512 A 
(Xi = 1-4), solving the implicit equations (^) and (||) results in H^^ k — 65.2 Gauss and 
fft^n = 77.0 Gauss. The ratio H^/H^ k = 1.18 and the difference - H%% k - 

11.8 Gauss. Table || presents other ferrofluids results: We note that the ratio £f*^™/-ff 
increases and the difference H^lt — Hcrit k decreases as Xi increases. It is better to use a 
ferrofluid with a low value of Xi to put experimentally in evidence the difference between 
H**$ and H%%K 



crit 

thick 
crit 
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6 Conclusion 



We get the dispersion relation of a layer of magnetic fluid of any thickness and viscos- 
ity, under a uniform vertical magnetic field. The critical magnetic field and the critical 
wavenumber of this instability are found to be thickness-dependent. The dispersion equa- 
tion, simplified into four asymptotic regimes, enables to explicit the expression of the 
critical magnetic induction of a thin film of ferrofluid. Near the onset of the instability, 
we show that the behaviour of the ferrofluid may be viscous or completely inertial; the 
behaviour which is manifested depends on the characteristics of the ferrofluid (contained 
in the parameter / proportional to v). In order to put in evidence experimentally the 
fact that the critical magnetic induction depends on the thickness, it is better to use a 
ferrofluid with a low value of the initial susceptiblity Xi- 
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Figure 1: Curve of marginal stability; thick-inertial layer. 
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Figure 2: Ferrofluid layer in the vacuum and in a magnetic induction. 
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Figure 4: Viscous and inertial regimes, thick layer; / = (Z„/Z c ) 3//2 and e 
$-1. 
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Figure 5: Viscous and inertial regimes, thin layer; / = (l v /l t 
and a is the layer thickness. 
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Re > 1 


a/l x < 1 
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a/l x > 1 
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Table 1: The asymptotic regimes 



capillary length 


i c = ii/pg) 1 ' 2 


capillary time 


t c = (l c /g) 1/2 


capillary Laplace pressure 


Pc = 7 Ac 



Table 2: Capillary quantities 



ferrofluid 


/ 


H /Hcru k when begins the inertial regime 


EMG 507 


0.0068 


1.00002 


EMG 901 


0.039 


1.0008 


APG 512 A 


0.27 


1.04 


APG 314 


0.71 


1.25 


APG 067 


1.29 


1.8 



Table 3: Range of the viscous and inertial regimes in the thick case 



ferrofluid 


EMG 308 


APG 512 A 


EMG 900 


Xi 


0.3 


1.4 


4.2 


M sat (Gauss) 


600 


300 


900 


H%% k (Gauss) 


299.5 


65.2 


19.7 


(Gauss) 


318.2 


77.0 


25.5 


Tjthin 1 TTthick 
n crit 1 n crit 


1.06 


1.18 


1.29 


H%%-H%$* (Gauss) 


18.7 


11.8 


5.8 



Table 4: Critical characteristics of various ferrofluids 



20 



ferrofluid 


p (g/cm 3 ) 


7 (N/m) 


r] (mPa s) 


Xi 


EMG 507 


1.15 


~ 0.033 


2 


0.4 


EMG 900 


1.74 


0.025 


60 


4.2 


EMG 308 


1.05 


~ 0.04 


5 


0.3 


EMG 901 


1.53 


0.0295 


10 


3 


APG 512 A 


1.26 


0.035 


75 


1.4 


APG 314 


~ 1.2 


0.025 


150 


1.2 


APG 067 


1.32 


0.034 


350 


1.4 



Table 5: Physical data of ferrofluids (Ferrofluidics Corporation). 
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